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Abstract 



Hadron-hadron collisions at high energies are investigated in the Ultra- 
relativistic-Quantum-Molccular-Dynamics approach. This microscopic trans- 
port model describes the phenomenology of hadronic interactions at low and 
intermediate energies {^/s < 5 GeV) in terms of interactions between known 
hadrons and their resonances. At higher energies, ^/s > 5 GeV, the excitation 
of color strings and their subsequent fragmentation into hadrons dominates 
the multiple production of particles in the UrQMD model. The model shows a 
fair overall agreement with a large body of experimental h-h data over a wide 
range of h-h center-of-mass energies. Hadronic reaction data with higher pre- 
cision would be useful to support the use of the UrQMD model for relativistic 
heavy ion collisions. PACS: 24.10.Lx, 13.75.-n, 13.85.-t 
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I. MOTIVATION 



Relativistic heavy ion collision experiments at the BNL-AGS (Au(10.7 AGeV)+Au) and 
at the CERN-SPS (Pb(160 AGeV)+Pb) have yielded a large variety of fascinating data. 
Various observables like the strong J/\E' suppression, enhanced yield of intermediate mass 
dilepton pairs, enhanced (anti-)hyperon yields, the creation of antimatter clusters and strong 
transverse flow seem to indicate the formation of very dense and highly excited matter [|l|. 
Since these observables are connected in a non-trivial way it is a tempting task for theo- 
reticians to model high energy heavy ion collisions in a consistent way and simultaneously 
predict this wide range of observables from a few hundred MeV up to several thousand GeV 
per nucleon at LHC. 

Bear in mind, however, that up to now there is no unique theoretical description of 
the underlying hadron— hadron interactions, with their vastly different characteristics at 
different incident energies and in different kinematic intervals. Perturbative quantum chro- 
modynamics (pQCD) can be applied to describe hard processes, i.e. processes with large 
four-momentum, Q^, transfer. But pQCD is formally inappropriate for the description of the 
soft interactions because of the absence of the large Q^— scale. Therefore, low— collisions 
are described in terms of phenomenological models. 

Early on, multiple production of secondaries in relativistic hadronic collisions has been 
described within the hydrodynamic approach 0. Then Regge theory and multiperiph- 
eral models have been developed to understand the phenomenology of the soft interactions. 
They avoid the difficulties attributed to the statistical models. An inconvenient point of 
this approach is the large number of free parameters, which have to be fixed by compari- 
son to experiment. Subsequently, various QCD- motivated quark— parton models have been 
introduced. 

Consequently, a vast variety of models for hadronic- and nuclear collisions have been 
developed. They may be subdivided into macroscopic (statistical and hydrodynamical) 
models 0] and microscopic (string-, transport-, cascade-, etc.) models like, e.g. UrQMD [Q, 
which is applied in the present paper, FRITIOF 0, VENUS 0, QGSM [|], RQMD [| and 
others |pi]|-[TB[] including the parton cascade approach . In the hydrodynamical (thermal) 
model one assumes local (global) equilibrium - the dynamics is characterized by the equation 
of state employed. The microscopic models describe subsequent individual hadron— hadron 
collisions. 

For low and intermediate energies hadron— hadron and nucleus— nucleus collisions are 
described in terms of the interactions between hadrons and their excited states, resonances, 
i.e. on the (quasi-)particle level. At high energies the quark and gluon degrees of free- 
dom cannot be neglected. Then the concept of color string excitations is introduced with 
their subsequent fragmentation into hadrons. In lead— on— lead collisions at the full SPS 
energy one finds in the UrQMD model that the ten most frequent hadron— hadron collision 
types namely, A^vr, tttt, Att, A^A^, Trp, A^A, vrif, ttt], ttu, Kit (in decreasing order of frequency) 
describe only 50% of the total h-h collisions. The inclusion of an additional 120 h-h collision 
types allows modeling up to 90% of the collisions predicted in the UrQMD model, while 
several thousand different h-h combinations are needed to cover more than 99% of the total 
number of h-h collisions. 

Since only a few of these cross sections are measured, one relies heavily on extrapolations 
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(and transformation, e.g. via the detailed balance principle) of known processes. Therefore, 
here we want to present and analyze the detailed elementary h-h input used in the transport 
model UrQMD. The h-h predictions are a necessary basis for our understanding of the 
dynamics of the complex heavy ion reactions. 

This paper is structured as follows: a brief description of the basic principles of the 
UrQMD model is given in Sec. ||. Section |T| presents UrQMD-results of the different h-h 
cross-sections for different reactions and a comparison with the available experimental data. 
The additive quark model (AQM) is used to calculate unknown cross-sections. UrQMD is 
the first microscopic model which attempts to include the color coherent phenomena. The 
implications of the effects of color opacity and color transparency in the model are discussed. 
The treatment of the formation and decay of resonances and strings is described in detail in 
Sec. 0. The importance of finite size effects in the fragmentation of strings is demonstrated. 
Section discusses the generation of the transverse momentum of particles in the model. 
In Sec. |VI| several predictions of observables for elementary channels are presented, which 
are especially interesting for the upcoming proton-proton run of the NA49 collaboration. 
Finally, a summary and conclusions are given. 



II. THE URQMD APPROACH 

The UrQMD-model |^ is a microscopic transport theory based on the covariant propaga- 
tion of all hadrons on classical trajectories in combination with stochastic binary scatterings, 
color string formation and resonance decay. It represents a Monte Carlo solution of a large 
set of coupled partial integro-differential equations for the time evolution of the various 
phase space densities fi{x,p) of particle species i = N,A,A, etc., which non-relativistically 
assumes the Boltzmann form: 

dfi{x,p) _ dp dfi{x,p) ^ dx dfi{x,p) ^ dfi{x,p) . . . 

where x and p are the position and momentum of the particle, respectively, and St/j(x,p) 
denotes the collision (or rather source-) term of these particle species, which are connected 
to any other particle species fk- 

The exchange of electric and baryonic charge, strangeness and four momentum in the 
i-channel is considered for baryon-baryon (BB) collisions at low energies, while meson- 
baryon (MB) and meson-meson (MM) interactions are treated via the formation and decay 
of resonances, i.e. the s-channel reactions, t-channel reactions for MB and MM collisions are 
taken into account from ^/s > 3 GeV on increasing to the only MB, MM interaction type 
above ^/s = 6 GeV. For nucleus-nucleus collisions the soft binary and ternary interactions 
between nucleons can be described by the real part of the in-medium G-Matrix, which is 
approximated by a non-relativistic density-dependent Skyrme potential of the form 

^^'^ = E ^i^i - ^j) + E S{x, - Xj)5{xj - Xk) , (2) 

where Xa denotes the coordinate variable in the quantum phase space. The first term 
simulates the attractive potential of the NN-interaction, and the second one yields the 
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saturation. In addition, Yukawa and Coulomb potentials are implemented in the model. 
The potentials allow to calculate the equation of state of the interacting many body system, 
as long as it is dominated by nucleons. Note that these potential interactions are only used 
in the model for baryons/nucleons with relative momenta Ap of less than 2 GeV/c. For the 
hadronic collisions discussed here, the potential interactions are omitted. Further details of 
the application of the UrQMD model to heavy- ion reactions may be found in . 

This framework allows to bridge with one concise model the entire available range of 
energies from the SIS energy region (-y/i ^ 2 GeV) to the RHIC energy {y/s = 200 GeV). 
At the highest energies, a huge number of different particle species can be produced. The 
model should allow for subsequent rescatterings. The collision term in the UrQMD model 
includes more than fifty baryon species and five meson nonetts (45 mesons). The baryons 
and baryon resonances included in the UrQMD are listed in Table |. In addition, their an- 
tiparticles have been implemented using charge-conjugation to assure full baryon-antibaryon 
symmetry. Figures |l| and |^ depict the implemented meson multiplets: pseudo-scalar, vector, 
scalar, pseudo- vector and (not shown in the Figs.) the tensor mesons as well as the heavy 
vector meson resonances p(1450), p(1700), u;(1420), and ti;(1600). Extremely heavy meson 
resonances (m > 2 GeV) are not explicitly implemented, however they may be important 
when investigating, e.g. the dynamics of correlations in future experiments. 

All particles can be produced in hadron-hadron collisions and can interact further with 
each other. The different decay channels all nucleon-. A- and hyperon-resonances up to 
2.25GeV/c^ mass as well as the meson (e.g. K*) decays etc. are implemented. At higher 
energies we take advantage of the hadron universality and use a string model for the decay of 
intermediate states. The cross-sections of various hadronic processes as well as the formation 
and fragmentation of the strings are discussed in the subsequent chapters. 



A basic input into the microscopic transport models are the particle species and -energy 
dependent cross-sections of hadron-hadron interactions. The total cross-sections are inter- 



preted geometrically. A collision between two hadrons will occur if d < y (Xtot/Tr, where d and 
(Ttot are the impact parameter of the hadrons and the total cross-section of the two hadrons, 
respectively. In the UrQMD model the total cross-section cxtot depends on the isospins of 
colliding particles, their flavor and the cm. energy. However, partial cross-sections are 
then used to calculate the relative weights for the different channels. Only a small frac- 
tion of all possible hadronic cross-sections has been measured. In the following sections, we 
compare the UrQMD cross-sections with experimental data. If no data are available, the 
additive quark model and detailed balance arguments are used to extrapolate such unknown 
observables. 



III. CROSS-SECTIONS 




A. Baryon-Baryon Cross-Sections 



The total BB cross-section of the reaction A + C ^ D + E has the general form 



aZfiV^s) oc i2So + + 1)^ ^ \M\' 




(3) 
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with the spins of the particles, Si, momenta of the pairs of particles, < pij >, in the two- 
particle rest frame, and the matrix element The matrix element |A^P, however, can 
take on a very complicated form and may be in general a function of of all the particle's 
quantum numbers as well as it's momenta and the cm. energy. 

If high quality experimental data on the respective cross section exists, a phenomeno- 
logical fit to the respective data is by far the most accurate approach for implementing the 
cross section. Otherwise, we have to rely on simplified assumption for the matrix element 
and employ general symmetries, like the principle of detailed balance (see section |IV A| ) . 

Let us start by investigating the total cross-section of proton-proton collisions from a 
beam momentum of 0.1 GeV/c up to 10"^ GeV/c as shown in Fig. ^ The total and inelastic 
cross-sections of the pp reaction are well measured in this energy region [jl5|. One finds 
a complex structure in this cross-section: local minima at 700MeV/c {Ecm ~ 2 GeV) and 
100 GeV/c {Ecm ~ 10 GeV), the maximum at 2 GeV/c {Ecm ~ 2 — 3 GeV) and a rise above 
lOOGeV/c. Note that the steep rise in the data below 300MeV/c is due to soft Coulomb 
interaction of the protons and is taken care of via the potential interaction. 

The structure in the pp cross-section is mainly due to the inelastic channels which are 
shown in Fig. ^. One clearly sees the A excitation with its increasing cross-section at low 
energies. The different partial cross sections depicted in Fig. ^ are discussed in section |1V A . 

A detailed comparison of prominent outgoing channels is depicted in Fig. ^ Here 
we show calculations of exclusive {pp —>■ mpp) and inclusive {pp mX) cross sections 
for the production of neutral mesons m = 7r^,ri,p^,uj as a function of the excess energies 
e = ^/s — ^/sih■ Here ^/s^ is the energy of the production threshold calculated as Y^ith = 
2mp + nim with the proton mass rUp and the pole mass of the meson rUm- In the case of 

mesons we count only those with masses within ±100 MeV around the pole mass to 
compare with data. The exclusive rj production just above threshold [|TB[ is overestimated 



by a factor of about two. Note that above 3.5 GeV the exclusive cross sections become 
less important because the string and multiple decay channels open and allow for multiple 
resonance production. For upcoming GSI-SIS experiments relevant y/s values are below 
4 GeV, where so far no data on the inclusive channels are available and therefore rely 
heavily on extrapolations. 

In Fig. 1^ the cross section oi pp ^ p + X reactions is shown. Good agreement with the 
data is found over a large energy range. 

At higher energies, the contributions of the different nucleonic resonances decrease and 
give way to the excitation of color strings, which is the dominant process at high energies in 
our model. The total cross section above the resonance region is given by the CERN-HERA 
parameterization [|T^ as shown in Table 

Partonic pQCD scattering is not included into the UrQMD model in the present version. 
The difference between the total and the elastic cross-section is taken as the inelastic cross- 
section. 



B. Meson-Baryon Cross-Sections 

The MB cross-sections are dominated by the formation of s-channel resonances, i.e. the 
formation of a transient state of mass m = y/Shhi containing the total cm. energy of the 
two incoming hadrons. On the quark level such a process implies that a quark from the 
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baryon annihilates an antiquark from the incoming meson. Below 2.2 GeV cm. energy 
intermediate resonance states get excited. The total cross-section of these reactions are 
given by the expression: 



^(v^) = {jB,mBjM,mM\\JR,MR) 

R=A,N* 

^ R^MB^tot 



{2Sb + 1){2Sm + 1) 



'>L (M^.-v^)^ + rL/4 ' 

which depends on the total decay width Ftot, on the partial decay width Tji^mb and on 
the cm. energy y/s. At higher energies the quark-antiquark annihilation processes become 
less important. There, t-channel excitations of the hadrons dominate, where the exchange 
of mesons and Pomeron exchange determines the total cross-section of the MB interaction 

Figures |^ and ^ show the cross-section of pion-proton reactions at different energies. 
In Fig. (vr^ + p) one probes predominantly the creation of the A"*""*" (A*"'"+) resonance. 
Note that the low energy s-wave rep scattering is not included into the UrQMD fit. The 
resonance peak ai p = 1.5 GeV/c is from the A(1900 — 1950) resonances. In comparison 
Fig. 1^ (vr" + p) depicts many strong uncharged non-strange baryon resonances, e.g. the 
A°(1232), A°*(1620), . . . , A^*(1535), etc. The total cross-section in the intermediate energy 
regime is therefore the sum of the individual excitation modes of baryon resonances, the 
s-wave at lower energies is left out. 

Let us now investigate collisions of strange mesons with baryons. For qs mesons strange 
s-channel resonances can be formed on non-strange baryons due to the annihilation of the 
g-quark. A comparison of these processes from UrQMD with the experimental data is shown 
in Fig. 1^. The formation of hyperon resonances is clearly visible at lower energies, while 



the universal t-channel reaction dominates the high energy tail. Fig. |I0| shows the cross- 
section of i^^-mesons {us) on protons. In this case, the formation of resonances is forbidden, 
since the s-quark cannot be annihilated by non-strange baryons (For strange baryons the 
formation of resonances is still possible). Here we use only the elastic channel and the t- 
channel excitation of both particles. The cross section at very high energies is given by the 
CERN-HERA parametrization as shown in Table y. 



C. Meson-Meson Cross Sections 

Due to the fact that the experimental preparation of meson beams and targets is re- 
stricted to tt's and K's, only very little is known about MM collisions in general. For the 
description of heavy ion collisions the importance of this channel increases with energy: at 
1 A GeV beam energy we find that the production of new hadrons (mostly pions) is only 
a ten percent effect. At AGS energies (10 A GeV) the amount of mesons roughly equals 
the number of incoming nucleons. Going on to the SPS (160 A GeV) the picture changes 
drastically: The produced particles dominate the reactions, while the incoming nucleons 
have dropped to a 15% admixture in particle density and multiplicity . 



To describe the total meson-meson reaction cross-sections, we make use of the additive 
quark model (see below) and the principle of detailed balance, which assumes the reversibility 
of the particle interactions. 
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Fig. |TT| compares the calculated cross-section of ti^ti~ scattering to experimental data 
T9| . The spectrum is dominated by the formation of the p with a mass of 770 MeV, the 
other two small peaks belong to the /o(970) and /2(1270) resonances. The /o(970) resonance 
is not visible in the data since the experimental analysis of meson-meson scattering is model 
dependent. 

Figure [l^ shows the implemented elastic tt^ti^ cross-section, which remains constant in 



the whole energy range where data are available pn| , pT |. 

Strangeness production in the meson-meson channel is possible, e.g. via the reaction 
TT+TT^ KK as shown in Fig. 113. 



Rescattering of strange mesons is implemented via the resonance formation (cf. Fig. ^ 
dominated by the K* resonance l2^), or elastically as depicted in Fig. |l^ |l23[| . 



Finally, we predict the cross-sections of tt+vt", t^^ and 7r"'"?7-reactions (Fig. |I6|) which 
are of utmost importance for the production of thermal photons and dileptons. At higher 
energies other meson resonances can be formed. To model MM interactions above the 
resonance region ^/s > 1.7 GeV we use the rescaled total vrp cross section: 

^tT(v^>l-7GeV) = aS(v/^)^ . (5) 

'^AQM 

This is justified, since at high energies the total cross section is given by quark counting. In 
the energy region from ^/s > 1.7 GeV to y/s < 6 GeV s-channel interactions are taken into 
account, while from ^/s > 3 GeV on t-channel excitation of both mesons becomes the MM 
interaction process of increasing importance in the model. 

The cross section for high energetic reactions are taken from the AQM-rescaled n^p 
paramatrization by the CERN-HERA group (see Table |ID. 



D. Antibaryon— Baryon Cross-Sections 

The physics of baryon— antibaryon interactions has been an area of much theoretical 
and experimental activity for a rather long period. It is well-known that at energies 
Piab < lOOGeV/c an important contribution to the total interaction cross-section comes 
from the process of annihilation, where only mesons are left in the final state. Though 
the earlier experiments on pp— annihilation revealed a number of differences from the non- 
annihilation channels, it is not clearly understood whether these differences arose simply 
due to the kinematic restrictions on the available phase space, or whether they are related 
to dynamical differences between the non-annihilation and annihilation mechanisms. The 



experimental results obtained in by comparison of pp with non-annihilation pp inter- 
actions at 32 GeV/c support the conclusion of equivalence of pp and non-annihilation pp 
interaction processes. 

Still, the nature of the baryon annihilation is subject to theoretical discussions. In 
the framework of the quark models based on topological 1/A^ expansion, the annihilation 
process is associated with the annihilation of string junctions, i.e. the point where strings 
are connected, such that the baryons have a y-form. In this case three qq strings are formed. 
The theory also allows for diagrams where the string junctions and one or two of the valence 
quarks can annihilate, corresponding to the formation of two strings or one string. Other 
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theories consider an annihilation mechanism without invoking the string junction hypothesis. 
This intriguing question has yet to be clarified (for review see, e.g. [25|, and references 
herein) . 

To avoid the difficulties attributed to these theoretical approaches, UrQMD is adjusted 
to known experimental data. The total pp cross-section is shown in Fig. as well as the 
annihilation and the elastic cross-sections. The UrQMD parameterizations depicted by lines 
are taken from Koch and Dover [EH]: 



a: 



pp 



S 



^0 — 



{s - so)2 + A^so 



+ B 



(6) 



120 mb. So 



4m%, A 



50 MeV and B = 0.6. The np cross-section does not 



with aQ 

differ significantly from the pp cross-section ||2^, hence they are set equal in the UrQMD 
model. 



At higher energies, CERN-HERA parameterizations |T5[ are used for the total and elastic 
channel: 



^ZAp) = A + Bp" + Cln2(p) + Dln(p) 



(7) 



with the laboratory momentum p in GeV/c and the cross-section a in mb. The parameters 
of the fit are listed in Table 

Below piab < 5 GeV/c the following parameterization is used: 



Crtot(p) 



75.0 + 43.1p-i + 2.6j9-2 
271.6 exp(-l.lp2) 



3.9p 



0.3 < j9 < 5 
p < 0.3 



(8) 



) 31.6 + 18.3J9-1- l.lp-2_3.8p . 0.3<p<5 
^^^^P^ = \78.Q ; p<0.3 

The sum of annihilation and elastic cross-sections do not yield the total cross-section: 

Aa = (Jtot - CTel - CTann (lO) 

In UrQMD this difference Aa is interpreted as the diffractive cross-section which describes 
the excitation at least one of the collision particles to a resonance or to a string via Pomeron 
exchange. 

The annihilation of baryon— antibaryon pairs proceeds in the UrQMD model accord- 
ing to rearrangement diagrams. Here the formation of two gg-strings of equal energies in 
the cm. system is assumed while the remaining constituent quarks are rearranged into 
newly produced hadrons. The generalization of the pp cross-section towards all possible 
antibaryon— baryon collisions can be done in different ways: 

1. The anti-baryon baryon cross-section at a given cm. energy ^/s is equal to the pp 
annihilation cross-section at the same ^/s^. 

(^Bsly^ = (^pp\^ (11) 
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2. The anti-baryon baryon cross-section at a given relative momentum prei is equal to 
the pp annihilation cross-section at the same Prei (Fig. 



Prei 



a. 



PPlPrel 



(12) 



In the UrQMD we have chosen the first parameterization. Since (Xann oc the 
annihilation cross-section drops rapidly with rising particle mass. The different treatment 
this cross-section can lead to systematic shifts in the antibaryon distributions for massive 
systems. Therefore, anti-baryon production off nuclei may be used to solve this ambiguity. 



E. The Additive Quark Model (AQM) 

Unknown cross-sections are calculated on the basis of the Additive Quark Model (AQM) 
which assumes the existence of dressed valence quarks, interacting very weakly 



inside of the hadron. At the phenomenological level, the AQM gives a correct quantita- 
tive and qualitative description of, e.g., the asymmetry of cm. spectra of secondaries in 
meson— nucleon and photon— nucleon reactions. The AQM has predicted also the important 
role of resonances for the multiple production of particles in hadronic interactions. To apply 
the AQM to calculations of heavy ion reactions, one needs to know the cross-sections of the 
quark interactions, which can be evaluated from the quark masses. Then, the unknown total 
cross-section of the high energy reaction can be calculated assuming a 40% reduced s-quark 
cross-section (compared to that of u- and (i-quark). The elastic cross-section is derived from 
Regge theory [pO|]: 

3 

c^eiastic = 0.039al^^^ [mb] , (13) 

where 

atotai = 40 (2\™^+™^ A _ 0.4-^) (l - OA-^) [mb] , (14) 



3/ V 3 — nil J V 3 — m2 

where rrii = 1(0) for particle i being a meson (a baryon) and is the number of strange 
quarks in particle i. This formula results from the high energy reactions, therefore there is no 
difference between antiparticles and particles. For BB-reactions no additional energy depen- 
dence is employed in collisions involving strange baryons. Non-strange baryon cross-sections 
are not treated via the Additive Quark Model, they have an explicit energy dependence in 
line with experimental data. The MB and MM cross-sections are rescaled via: 

/ cr7rAf(v^) AQM /1r-^ 
(^X^X2WS) = AQM~^X,X2 ■ (15) 



■N 



The hyperon-nucleon cross-section, which is taken in the UrQMD model from the Additive 
Quark Model, is in good agreement with the data above piab = 300 MeV (cf. Fig. |1^). The 



total cross-sections calculated for baryon-baryon, meson-baryon and meson-meson interac- 
tions are listed in Tables [m|— |V|. 
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F. Color Fluctuations, Color Opacity and Color Transparency 



Quantum chromodynamics has important applications of the dynamical role of color 
degrees of freedom to the strong interactions at ultrarelativistic energies (for a review, see 
|31| , |32[| and references therein). The theory is presented in Ref. ^1|] in detail, here we 



just sketch the main ideas of color optics and coherent phenomena in high energy physics. 
Hadrons are dynamical objects which come in Fock space of configurations of very differ- 
ent spatial sizes. At high energies, incident hadronic quark— gluon configurations can be 
considered frozen as a result of Lorentz time dilation. Due to the long coherence length at 
such high energies one can apply geometrical color optics. Small objects produced in hard 
processes with high momentum transfer then have reduced interaction cross sections. In 
processes with moderate such compact objects, which are a coherent superposition of 
eigenstates of the QCD Hamiltonian, should gain size. When the quark— gluon configuration 
is large, it will lead to an increased interaction cross-section of the hadron with the nuclear 
medium. 

Therefore, the fluctuations of the hadron's spatial extent give rise to the color trans- 
parency and color opacity phenomena: When a small object is produced, it interacts only 
very weakly with other hadrons due to color screening. Moreover, since - at sufficiently high 
energies - the small-sized configuration of this object is frozen, the nuclear medium appears 
to be transparent for such hadrons (color transparency). In contrast, hadronic configura- 
tions which are larger than average interact with larger cross-section, giving rise to color 
opacity. 

Nucleus-nucleus collisions provide a tool to investigate the effect of color transparency, for 
instance, in the production of leading nucleons. On the other hand, the complementary color 
opacity effect, i.e. large-sized configurations, can cause stronger stopping and significant 
fluctuations in the transverse energy of secondaries in central reactions. 

A first step to investigate these QCD effects within a microscopic transport model is 
made by incorporating the color fluctuations in the elementary hadron— hadron reactions in 
the UrQMD model. Thus, one needs to know the probability -P(o") that a given configuration 
interacts with a nucleon with a total cross-section o. It is convenient to consider moments 
of the distribution: 



(aO) = y daa°P(a) = 1 , (16) 
(a^) = f daa^P{a) =a ,etc. , (17) 



where a denotes the average cross-section. The second moment < > can be determined 
from the diffractive dissociation experiments. In addition, further information can be ob- 
tained from QCD, which implies [ pi]] : 

P((t) oc , (18) 

for (T — *■ 0, where A^^^ is the number of valence quarks. Thus, for the nucleon and the pion 
distributions it follows for a ^ 0: 

P^(a)oca , (19) 
PTr{(y) oc constant . (20) 
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Prom these arguments, P{cr) can be construct. Fig. ^ shows the resulting broad -P(cr) 
distribution for proton projectiles and the even broader one for the pions [0. 

The effect of the color fluctuations on proton-proton interactions at different impact 
parameters, b, is shown in Fig. |2l| for the UrQMD model. The charged pion multiplicity 
distribution decreases monotonically with rising b. A^7r(&) has a non- vanishing tail for b > 
1.1 fm, in contrast to the abrupt geometrical edge of the distribution as calculated in the 
static geometric models without color fluctuations. 



IV. THE REACTION CHANNELS 



A. Resonances 



The production and decay of resonances is the most important h-h reaction channel 
below ^/s = 5 GeV for BB and 3 GeV for MM and MB reactions. Baryon resonances are 
produced in two different ways, namely 

i) hard production: N+N^ AN, AA, N*N, etc. 

ii) soft production: vr^+p^ A*^, K^+p^ A*... 

The formation of s-channel resonances is fitted to measured data, e.g. in the reaction 
A + C ^ D + E we use the general form 

af,f(v^)oc(25z, + l)(2^i, + l)^^i|A<(m^,m^)P , (21) 

\Pa,c) s 

with the spins of the particles. Si, momenta of the pairs of particles, < Pij >, in the two- 
particle rest frame, and the matrix element |A^(m£), m^;)^, which here depends only on the 
masses of the outgoing hadrons, rui. 

There are six channels of the excitation of non-strange resonances in the UrQMD model, 
namely A^A^ ^ A^Ai232, A^A^*, A^A*, A1232A1232, Ai232Ar*, and A1232A*. The A1232 is explic- 
itly listed, whereas higher excitations of the A resonance have been denoted as A*. For each 
of these 6 channels specific assumptions are made with respect to the form of the matrix 
element Ai, and the free parameters are adjusted to the available experimental data, when 
available: 



1. A^A^ A^Ai232 excitation: 



\M{V~s,m„m,)\^ = A ^ , (22) 



with TUA = 1.232 GeV, Fa = 0.115 GeV and A = 40000. Note that this form of the 



matrix element has been adjusted to fit the data shown in figure 
2. A^A^ ^ A^A^*, A^A^ ^ A^A*, A^A^ ^ Ai232Ar* and A^A^ ^ A1232A* excitation: 

1 



|A^(m3,m4)|^ 



A 



(m4 — 7713)^ (m4 + 7713) 



(23) 
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with A = 6.3 for A^A^ ^ A^A^*, A = 12 for A^A^ ^ A^A* and A = 3.5 for A^A^ ^ 
Ai232A^*- In the case of A^A^ A1232A* there are insufficient data available, therefore 
we use the same matrix element and parameters as in the case of A^A^ — >■ Ai232A^*- 
Since rris 7^ is valid for all above cases, the matrix element cannot diverge. 

3. A^A^ A A excitation: 

\M{ms,m,)f = A , (24) 



with A = 2.8. 



Figure |2^ shows the fit of the UrQMD pp — > A^Ai232 cross section to experimental 
measurements [Q. The measurements refer to the A^ + n exit channel and have been 
rescaled to match the full isospin-summed cross section. In the case of the exclusive A1232 
cross section the quality of the data and thus also the quality of the resulting fit is very good. 
For higher resonance excitations this is unfortunately no longer the case and additional 
measurements are needed to clarify the situation. One has to keep in mind, however, that 
the experimental extraction of exclusive resonance production cross sections is only possible 
via two- or three-particle correlations (e.g. a pion-nucleon correlation in the case of the A) 
which introduces large systematic errors, especially for broad resonances. 

In figure |2^ the UrQMD cross sections for the processes pp — >• PP1440, PP — PP1520 



pp PP1680 pp — > pPijoQ are compared to data One single parameter has been 



used to fit all four cross sections. The data situation is not as good as in the case of the 
A 1232 resonance, some ambiguities are visible which results in the quality of the fit being 
not as good as in the previous case. The parameters for the other classes are fitted in the 
same fashion. 

The cross section for exclusive A1232A1232 can be seen in figure The data points 
56| hint at a resonance like structure which cannot be reproduced with the UrQMD ansatz 



for resonance-excitation cross sections. However, the data deviates considerably from other 
cross sections for resonance excitation (e.g. A^A^ N*A*). Considering hadron-universality 
and the similarities in all other resonance excitation cross sections this casts a certain doubt 
on the accuracy of the measurement of the resonance-like peak. 

Figure |2^ shows the UrQMD fit for the exclusive A^g2oA^ production. The same matrix 
element is used for the entire class of A^A^ A*N reactions. In the case of exclusive A*N* 
production the matrix element has been fitted to the A^232-^r680 ^^^^ channel (see figure p6D . 
The extrapolation to the case of A*232^r52o production can be seen in figure For the 
exclusive pp AA* reaction class the data situation is unsatisfactory, therefore we used 
the same matrix-element as in the pp — >■ A^A* case. 

The decay of the resonances proceeds according to the branching ratios compiled by the 
Particle Data Group [0. The resonance decay products have isotropical distributions in 
the rest frame of the resonance. If the resonance decays into N > 2 particles, then the 
corresponding A^— body phase space is used to calculate their A^ momenta stochastically. It 
is necessary to note that a consistent description of angular momentum distributions points 
to a rather intricate problem of transport theory itself: If one considers the whole scattering 
interaction to be described by one single quantum mechanical process there are correlations 
between the final and the initial stage. For instance the angular distribution of the final 
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particles with respect to the axes of the incoming momenta in the CMS system. However, 
a fitting of the angular distributions to experimental data may conflict with the basic as- 
sumption of transport theories that the multiple scattering processes can be considered to 
be of Markovian type, i.e. after each scattering process or resonance formation the outgoing 
particles completely forget about their entrance channels. In the case of a spin resonance 
there is no preferred direction for the emission of the final particles, while for spin 1 (and 
other) the different magnetic quantum numbers are statistically occupied, so that also in 
these cases there is no preferred angle of emission^. 

All produced particles are able to rescatter within the nuclear medium, therefore the 
excitation of resonances by the annihilation of mesons on baryons included as depicted in 
Fig. for the reaction vr^ + p ^ A++*^*\ 

Also the 7T~ + p channel (Fig. ^ shows a rich structure of baryon resonances. The total 
meson-baryon cross section is given by formula There, the total and partial decay widths 
also define the inverse reaction, i.e. the different decay-channels of the respective resonance. 
Thus, the principle of detailed balance is applied. Based on this principle we calculate 
all resonance formation cross sections from the measured decay properties of the possible 
resonances up to cm. energies of 2.25 GeV/c^ for baryon resonance and 1.7 GeV/c^ in the 
case of MM and MB reactions. Above this energy collisions are modeled by the formation 
of s-channel string or, at higher energies (beginning at ^/s = 3 GeV), by one/two t-channel 
strings. In the strangeness channel elastic collisions are possible for those meson-baryon 
combinations which are not able to form a resonance, while the creation of t-channel strings 
is always possible at sufficiently large energies (c.f. Fig. ^ for the formation of hyperon 
resonances and Fig. |10| for the non-resonant channel). At high collision energies both cross 
section become equal due to quark counting rules. 

In more general terms, the principle of detailed balance can be derived by assuming time- 
reversal invariance of the interaction Hamiltonian and can be formulated in the following 
way: 

a{y x) pI Qy = a{x y) pi , (25) 

with p denoting the c.m.-momenta of the particles and g being the spin-isospin degeneracy 
factors. Thus, if the cross section of the reaction x ^ y is known, the back reaction y ^ x 
can be easily obtained. This principle is in UrQMD widely applied for the calculation of 
baryon-resonance absorption cross sections, such as A(1232) + ^ + A^. For a detailed 
discussion of the application of the principle of detailed balance to resonance absorption and 
i/s-dependent decay widths we refer to 

As was mentioned above, not only baryon-baryon and meson-baryon collisions have to be 
included in the proposed scheme. At high energies and in massive A A systems meson- meson 
collisions may dominate the multiple production of secondaries. Unfortunately, there are 
only few channels for which the experimental information exists, like the process of vr+vr" 



scattering (Fig. |Tl|), which is fairly described by the UrQMD model. 



^For a detailed discussion of the influence of non- markovian processes in the transport theory of 



heavy ion collisions we refer the reader to Refs. [p^, 34 1 . 
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B. Strings 



Both, for the high energy regime and for baryon— antibaryon annihilation we apply a 
string model (similar to the LUND model |3^) to describe the inelastic reactions. The 
constituents, quark and diquark (or anti-quarks), of the incoming hadron also define the 
predominant emission patterns of the events. The amount of stopping in nucleus-nucleus 
is strongly correlated to the detailed dynamics of the diquark in the hadronic medium [^. 
Recently different additional mechanisms of baryon number transport in nucleus-nucleus 
reactions have been investigated: Baryonic junctions as suggested by |^ which yield an 



enhanced hyperon and proton production cross section at central rapidities di-quark 
breaking due to interaction with the hadronic medium as predicted by The di-quark 



breaking component is also taken into account in the UrQMD model (di-quark breaking 
probability= 10%). However, the above mechanisms are of minor importance in the UrQMD 
approach since rescattering of the leading di-quark with hadrons is explicitly taken into 
account. 

Since gluons are massless particles with spin J = 1, the static strong interaction between 
quarks at small distances (r -C 1/Aqcd) may be described by a potential Vq oc —as/r. At 
sufficiently large distances the color field between two quarks or anti-quarks transforms into 
the color string. Denoting the string tension as k one defines the linear string potential 

Vi = k\zi-Z2\ , (26) 

between the quarks/di quarks located at zi and Z2, respectively. This form of the potential 
is chosen from heuristic considerations, based on the quark confinement hypothesis, and is 
supported by lattice QCD calculations pOj l. 



The transverse directions have not to be taken into account, since they are negligibly 
small compared to the longitudinal excitation of the hadron string. Hence, we get the 
dynamics of the quark system (with quark momenta pi and P2) from the Hamiltonian H 

H = \pi\ + \p2\ + k\zi - Z2\ , (27) 

which leads to the following equations of motion for the massless endpoints of the string: 

dPi dH 

^ = -^ = -sign(..-..). , (28) 

dzi dH . , , , „^ 

_ = +_ = s,g,>(p.) . (29) 

A change in momentum is directly related to the sign of {zi — z'j), while the direction of 
propagation is defined by the sign of the momentum pi of the quark. This results in a 
typical "yo-yo" type evolution of the quark system. 

If the momentum transfer is large enough, the excitation of the string may exceed some 
critical limit. After that it will be energetically favorable to break the string into pieces by 
producing gg-pairs from the vacuum. Each of the produced gg-pairs will have small relative 
momenta in their rest frame. Owing to the fact that the color string is uniformly stretched, 
the hadrons produced as a result of the string fragmentation will be uniformly distributed 
within the kinematically allowed interval between i/min = and ymax = In {s/m'^). 
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The probability of the pair production process has been calculated by Schwinger for the 
case of an infinite homogeneous electrical field. His result can be used to motivate the decay 
of QCD color field between the quarks. Note that QCD is a non-abelian theory, therefore 
the color field need not be homogeneous and it is definitely not infinite. This leads to the 
modifications of string decays which will be discussed later. 

The probability |Mp for the creation of a quark- antiquark pair with mass m in a color 
field with a string tension k is: 



|M| = constant x exp , (30) 



K 



where a typical value for k is 1 GeV/fm. This relation is motivated by Schwinger's QED- 
based result for particle-antiparticle creation in a strong electric field (see the discussion of 
formula ^ in section [IV B 1|) . The relative production probabilities of the different quark 
flavors are adjusted to e"'"e~-data: 

u : d : s : diquark = 1:1: 0.35 : 0.1 . (31) 

The production of charmed (and heavier) quarks is strongly suppressed in the string picture, 
hence they are exclusively produced in hard QCD processes. The elementary diquark is 
introduced to allow for baryon-antibaryon production in the string. A schematic view of the 
decaying string is shown in Fig. — a non-strange baryon decays into a hyperon, a kaon 
and a pion. 

To decide which type of hadron is produced from the quark configuration that is created 
in the color field we choose in the case of a produced 

(i) baryon — the octet and decuplet with equal probabilities; 

(ii) meson — the meson nonet with a probability proportional to the spin degeneracy 
and inverse mean mass m, 

25 + 1 

-P multiplet OC r . (32) 

\^/multiplet 

The singlet states are mixtures of uu, dd and ss. They are projected onto SU(3) hadrons 
with the flavor mixing angles from the quadratic Gell-Mann-Okubo mass formula For 
the scalar mesons this formula is not applicable, here an ideal mixing angle (tan(0) = l/\/2) 
is assumed (The mixing angles are depicted in Table [VI| ). 

The Field— Feynman fragmentation mechanism which allows the independent string 
decay from both ends of the string is used in the UrQMD model. The string break-up is 
treated iteratively: String — hadron -|- smaller string. The conservation laws are fulfilled. 
The diquark is permitted to convert into mesons via the breaking of the diquark link, thus 
transporting the baryon number into central rapidities. 

On both sides of the fragmenting string the new particles are formed randomly. If a res- 
onance is produced, its mass is determined by a Breit-Wigner mass distribution. The trans- 
verse momentum is assigned to this particle according to a thermal momentum distribution, 
resulting in a temperature of 170 MeV. After that the fragmentation function determines 
the fraction of the longitudinal momentum of the string transferred to the hadron. This 
procedure can be described in a covariant manner by the light cone variables defined as: 
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t ± z and = E ±p . 



(33) 



The light cone momentum Phadron given to the newly produced hadron is: 



Phadron ^fraction Ptotal 



The fragmentation of a baryonic string reads: 

P~ {qqqqq) = ^^actionP" (^J^ + (P" - %7actionP") ^ 

String Baryon String 



(34) 



(35) 



The main input is the fragmentation function which yields the probability distribution 
P (^fraction 5 ''^t)- This function regulates the fraction of energy and momentum given to the 
produced hadron in the stochastic fragmentation of the color string. For newly produced 
particles the Field- Feynman function \TL\: 



p{z, 



± 1 

fraction;' 



constant x (1 



^fraction / 



(36) 



is used. P{z) drops rapidly with increasing z (Fig. ^). Therefore, the longitudinal momenta 
of e.g. produced antibaryons (Fig. |3D|) and pions (Fig. ^T]) are small (they stick to central 
rapidities), in line with the experimental data. The rapidity spectra of these particles have 
a characteristic Gaussian-like shape, in contrast to the baryon spectra in pp, as it is clearly 



seen in Figure 30 



The proton is on average less stopped, since it is build up from the leading diquark in 
the string (leading particle effect). Fig. ^ compares the xp distribution of protons and A's 
for the Feynman scaling variable xp = 2p\\/ ^/s measured in pp reactions at 205 GeV/c. The 
data on leading baryons can only be reproduced when a modified fragmentation function is 
used for the leading baryons (cf. Fig. |29|, dashed curve). This leading baryon fragmentation 
function is of Gaussian form: 



P(4action) = constant X exp 



± 

fraction 



-by 



2a2 



(37) 



with parameters a = 0.275 and b = 0.42. 

It is obvious that modeling the momentum loss in elementary collisions has a strong 
influence on the rapidity spectra of particles produced in heavy-ion collisions. This can be 
seen if one compares the spectra of p, A's (created from the leading baryon) and A's or 
mesons (created from newly produced quarks) in lead-lead collisions at the SPS energies [Q. 



1. Finite Size Effects 

Let us now discuss finite size effects in the process of string fragmentation. The string 
is essentially a color field which connects two color charges, the [3] quark and the [3] di- 
quark or anti-quark, at the ends of the string. Quantum electrodynamics (QED) predicts 
spontaneous particle— antiparticle creation in strong electric fields |^2|,^. This effect should 
hold for particle creation in a strong color fields. Schwinger's QED result 
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[eEf ~ 1 

N=l 



1 \ 

m \ 



47]-3 



-N-K — I 



(38) 



is often adopted to the case of color fields by equating \eE\ with the string tension /t. 

However, for strings, several important assumptions, which lead to Schwinger's result 
are not fulfilled: Firstly, the color field is not infinitely extended. It is bound radially by the 
interaction length of the gluons and longitudinally roughly by the [3] and [3] endpoints. 
Secondly, the two endpoints of the string move (with close to the speed of light) in opposite 
directions. Finally, the assumption of a constant field strength k seems to be fulfilled along 



the string, but may become invalid in hot and dense matter K 



The influence of the finite radial size of the string can be studied by solving the Dirac 
equation for the (newly produced) pairs in a finite volume. The field is restricted to a 
cylindrical volume of length L and and area vri?^, where R is the cylinder radius. Along 
the cylinder axis a homogeneous color field is assumed. The boundary condition on the 
surface of the cylinder leads to enclosure of the color charges. Taking L —>■ oo neglects the 
longitudinal direction. 

The constraint on the cylinder surface is given by linear MIT-boundary conditions |^ 
thus we have to solve the equations: 








(x G dV) 



(39) 
(40) 



The MIT model allows for an analytic solution |4^. The occupation numbers in the limits 
t — > ±cx) yield the pair creation probability as |^ 



R 



{eEf 

47]-3 



E 

N=l 



exp 



-Nir 



m 
eE 



\ 2t^N ^ 
X < „ y ^ exp 

n,fj, 



\eER'^ 



-Nir 



eER^ 



(41) 



Thus the string radius is directly related to the pair creation rate Wr, with /c„^ being the 
j^th momentum eigenvector in the solution of Eqs. (^ and ^0|) for a given projection /i of 
the corresponding Bessel functions. 

The expression in curly brackets gives the deviation from the Schwinger formula. For 
i? — > oo this second factor vanishes and one is back to the infinite case. This result allows 
to calculate the strangeness suppression fs{R) = Wr{s) /Wr{u) for different string radii as 



shown in Fig. 33 



The overall pair creation rate has been calculated as a function of the longitudinal size of 



the color field for the corresponding Dirac equation in terms of confluent hypergeometric 
functions. The resulting particle production rate per volume dV and time interval dt, 
dN / {dV dtdpx) at = is shown in Fig. ^ as a function of the distance from the 
string center, for different string lengths L. A clear depletion of particle production near the 
endpoints of the string endings is visible (Eq.|38|). For short strings the particle production 
gets enhanced toward the string center. 

The above discussed finite size effects have a strong influence on the results, especially 
on the production of heavy quarks and particles like anti-baryons, -hyperons, etc. Up to 
now these corrections are included only in the non-strange antibaryon sector. This effect 
is very important especially in the case of the A-particle, since it is mostly created at the 
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string ends in the fragmentation of a leading proton. Indeed, it has been reported that most 
of the transport models which use a string fragmentation scheme based on the Schwinger 
formalism tend to overestimate the A's [HHI. Therefore, this finite size effects shall be included 



for strange baryons, in order to correct the overestimated number of A's and A's. 



V. THE GENERATION OF TRANSVERSE MOMENTUM 

In hot and dense nuclear matter, most hadrons suffer interactions from the many sur- 
rounding particles. As a result, the effective mass may change with density. Many dynamical 
properties of hadrons are modified in the medium. In-medium two-body scattering cross- 
sections may differ from the free space values. Those effects can be studied in the framework 
of relativistic transport theory, i.e. the relativistic Boltzmann-Uehling-Uhlenbeck (RBUU) 
equation. This type of transport equation has been used extensively to the study of rela- 
tivistic heavy-ion collisions and turned out to be very successful. The UrQMD approach 
uses an analytical expression for the differential cross-section of in-medium A^A^ elastic scat- 
tering derived from the collision term of the RBUU equation to determine the scattering 
angles between the outgoing particles in elementary hadron-hadron collisions. It is assumed 
that the angular distributions for all relevant two-body processes are similar modified in 
an analogous manner. They are approximated by the differential in-medium A^A^ elastic 
scattering cross-section: 

aNN^NN{s,t) = —\^[D{s,t) + E{s,t) + (s,t < — >u)], (42) 
[Zn )^s 

with the direct term 

Dis, t) = ^^^,it - Am*y + (f^(2^^ + 2st + - 8m*^s + 8m*-) 

and the exchange term 

, / r *2\ ^{SnNJ I a *2 j.\ *4j. 

X (s — 6m ) — — 5T(4m — s — t)m t 



{t — m1){u — ml) 



+ 



^ ^ ^ ?,m*'^{Am*^ - s-t){Am*'^ -t) ^ 3t{t + s)m*^ 



2(t — m^)(u — m^) 2(t — m'^)(u — m^.) 

- 4m*^s - 10m*H + 24m*^ {t + s)^ - 2m*^s + 2m*H 



4(t — ml) [u — m^) 4(t — m^) {u — m^) 

^ ^^, ?>m*\t + s-Am*^){t + s-2m*^) , 3m*\e-2m*H) ^ 
+ y i^t - ml){u - ml) + (t - ml){u - m^)^' ^^'^ 



The (pseudo-)scalar and vector coupling constants are g]^^ = 1.434, g^^ = 6.9, and 
7.54 and m* is the in-medium m 
single-particle energy is given by 



7.54 and m* is the in-medium mass, s, t, u are the Mandelstam variables. The in-medium 
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= ^p*2 + . (45) 

The formula for the differential cross section of in-medium NN elastic scattering can be 
used for elementary hadron-hadron collision if it is scaled by 



s ^ s-{ml + mlf + 4m*' , (46) 

where m\ and denote the effective masses of the incoming hadrons. Furthermore, we 
take into account finite size effects of the hadrons and part of the short range correlation by 
introducing a phenomenological form factor at each vertex. For the baryon-baryon-meson 
vertex the common form 

Fbbm = .2 + ■ (4"^) 

is used, where A m is the cut-off mass for meson M. 

The total energy and the masses of the incoming hadrons serve as input for calculating 
the angular distribution. It is worth to stress again that Eq. ( ^21) is used to calculate only the 
angular distributions for all elementary elastic two-body processes but not the corresponding 
total cross sections. The inverse slopes ('temperatures') as calculated in the UrQMD from 
the transverse momentum spectra of pions - by fitting the l/mtdN/drrit distribution with an 
exponential - compare well with values extracted from thermal fits to dataQ |^0| (Fig. |35|) . 



VI. PARTICLE YIELDS, LONGITUDINAL AND TRANSVERSE SPECTRA 

The UrQMD model reproduces nicely the total, elastic and inelastic cross-sections of 
numerous hadronic reactions. The model also predicts the particle multiplicities (i.e. the 
inclusive cross-sections) as well as the (Lorentz- invariant) cross-sections, which may come 
in the form of xp-, rapidity-, or transverse momentum distributions. The abundances of 
the most important particle species produced in pp collisions at 12 GeV/c are listed in 
Table [VTTI - the model predictions are in line with the data within 15%. The yields of various 
particle species produced in pp collisions at 205 GeV/c {^/s ^ 19.7 GeV) [^] are listed in 
Table |V111| . It is easy to see that the model predictions lie generally well within the 10% of 
the data except for the strange baryons. Table presents the results on particle production 
in pp interactions at ^/s = 27 GeV [Q. Again, the agreement between the experimental and 
theoretical results is good. As discussed above, the production of A(A)'s is overestimated 
by a factor of 2-3 due to the neglect of finite size effects in the strings. 

The rapidity spectra of pions and Xir-distributions of baryons, as obtained from the 
UrQMD model for pp interactions at 205 GeV/c (Fig. |30| and respectively) have been 
discussed already in the previous section. Correlations between the transverse and longitu- 



dinal momenta of charged pions, produced in the same reaction, can be studied |j52|. The 



^Temperatures have been extracted from a statistical model fit to the particle yields in pp, pp and 
e+e" . 
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transverse momenta of both positively and negatively charged pions have evident minima at 
xp = SiS shown in Fig. ^ (the "sea-gull" effect). Then the transverse momentum increase 
nearly linearly with the rising longitudinal momenta. The agreement with the experimental 
data is fair. The single event correlation between the transverse momenta of the produced 
pions and the multiplicity of the negatively charged hadrons in the event is presented in 
Fig. |3^. The measured distribution is reproduced nicely for tt+'s , while for the 7r~'s, 20% 
deviations of the calculated spectrum from the data are evident at high [uh- > 5) event 
multiplicities. 

To probe the ability of the UrQMD model to describe hadron-hadron collision even at the 
today highest available bombarding energies for nuclei, we compare the calculated He+He 
collision at ISR with data |^ as shown in Fig. |3^. 

It is not surprising, that the light helium system is transparent. A baryon free area of 3 
units in rapidity is produced. The UrQMD model prediction describes the data fairly well. 
The calculated particle yields can be increased by 15% if one includes multiple jet creation 
into the model description . 



The UrQMD model reproduces the cross-sections and spectra of particles in hadronic 
collisions fairly well. Since hadronic interactions build up the basic input to simulate the 
hadron-nucleus and nucleus-nucleus interactions in the model, it is interesting to examine 
the applicability of the UrQMD model to these reactions. The full comparison with the 
experimental hA and AA data is an ongoing program, which is not completed yet . 



The ability of the model to reproduce e.g. dilepton yields in pp collisions, which is 
of interest in high energy physics is shown here as an example. Fig. ^ shows the UrQMD 
calculations of the dilepton spectrum for p+Be (which serves as a substitute for pp reactions) 
at 450 GeV/c. Dilepton sources considered here are Dalitz decays (vr", rj and u) and direct 
vector meson decays (p, u and 0). Dalitz decays of heavier meson and baryon resonances are 
included explicitly via their emission of p mesons (assuming vector meson dominance). To 
avoid double counting, the p mesons from r/'s, and cj's are excluded from the p contribution. 
Pion annihilation is included dynamically into the contribution of decaying p mesons via the 
channel tc^tt^ — > p — »■ e'^e~ . The calculation of dilepton yields without modifications of the 
p mass pole is compatible with the CERES data |5H . 



VII. SUMMARY 



The basic hadronic interaction processes incorporated into the UrQMD approach are 
given. The implemented cross-sections of various hadron-hadron collisions, as well as their 
extrapolations into the high energy region are presented. The model treatment of excitation 
and decay of intermediate objects like resonances and strings is described. The importance 
of the finite size effects for the process of string fragmentation is discussed. The UrQMD 
model is a microscopic transport model which allows to study optionally color fluctuations, 
i.e. color coherence phenomena, as Color Transparency and Color Opacity, as well as the 
expansion of small wave packets from the point of the production. The model predictions are 
compared with the available experimental data on particle yields in hadronic interactions 
for a wide range of cm. energies. The UrQMD model treats the elementary processes 
reasonably well. More accurate data on proton-proton, proton-neutron, as well as meson- 
baryon collisions are needed to improve the extrapolation to nucleus-nucleus interactions at 
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high energies, where already in the present model - about 10^ elementary hadron-hadron 
reactions are possible. 
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FIG. 2. Implemented mesons: scalar mesons (left plot) and pseudo- vector mesons (right plot). 
fi and /( are the states /i(1285) and /i(1420), respectively. 
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FIG. 3. The total cross-section of collisions vs. the laboratory momentum piab of the incident 
particle. Data are taken from p5[|. 



27 



UrQMD l.a 




inelastic 

p p — iN A 

p p — *N A 

— pp^n'a 
H — P P — 

— p p — >N string 



(GeV) 



FIG. 4. The inelastic cross-section of pp collisions vs. the laboratory momentum piah and 
cross-sections of the various inelastic channels. 
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FIG. 5. Cross section for the production of neutral mesons in pp. The inclusive and exclusive 
meson production is compared to data by |16| 
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FIG. 6. Cross section for the production of anti-protons in pp as a function of cm. energy. The 



UrQMD calculation is compared to data by |57| 
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FIG. 7. The total cross-section of n^p interaction vs. laboratory momentum piab- Data are 
taken from [^. Note that the low energy s-wave n^p scattering is not included into the UrQMD 
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FIG. 8. The same as Fig. |^ but for vr p interaction. Data are taken from [15|. Note that the 
low energy s-wave ■7T~p scattering is not included into the UrQMD fit. 
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FIG. 9. The same as Fig. but for K p reaction. Data are taken from [15| 
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FIG. 10. The same as Fig. but for K~^p reaction. Data are taken from 
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FIG. 12. The same as Fig. but for TT^vr^ scattering. Data are taken from ||2^ (open squares) 
and from E^ (open circles). 
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FIG. 15. The same as Fig. 14 but for tt K reaction. Data (open squares) are taken from 
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FIG. 17. The pp cross-section as compared to the experimental data on total (open circles), 
elastic (open squares), and annihilation (open triangles) cross-sections. Data are taken from p5| . 
The diffractive cross-section is assumed to be a difference between the total cross-section and the 
sum of the elastic and annihilation cross-section. 
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FIG. 18. Extrapolation of the pp cross-section towards unknown antibaryon— baryon reactions. 
Here the cross-section of the antibaryon— baryon interaction is equal to the pp cross-section at the 
same relative momentum. In the UrQMD model we take the antibaryon— baryon cross section 
equal to the pp cross-section at the same center-of-mass energy. 
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FIG. 19. Total A-p cross-section vs. laboratory momentum of the A. The UrQMD results are 
given by the Additive Quark Model. Data are taken from |15|. There seems to be indication for a 
resonance at ^/s = 2.1 GeV. This could be a 6q molecule or a di-baryon with s = 1. 
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FIG. 20. Probability distribution of nucleon (dashed curve) and pion (solid curve) cross-sections 
as predicted by |3l||. 
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FIG. 21. Charged pion multiplicity in pp collisions at ^/s = 27 GeV with (dotted-) and without 
(solid-line histogram) color fluctuations for different impact parameters. 
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FIG. 22. UrQMD fit for tlie exclusive A1232 production in proton-proton reactions compared 
to data |56|. 
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FIG. 23. UrQMD parameterization for exclusive pp* cross sections. Only one parameter was 
used to describe all available cross section data [^] . 
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FIG. 24. UrQMD fit for the exclusive A1232A1232 production in proton-proton reactions 
compared to data ||56[. 
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FIG. 25. UrQMD fit for the exclusive Ai92oA^ production in proton-proton reactions compared 
to data |56|. 
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FIG. 26. UrQMD fit for tlie exclusive Ai232iV*goQ production in proton-proton reactions 



compared to data |56|. Tlie matrix element for all Ai232.^x reactions is extracted from this fit. 
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FIG. 27. Comparison between the UrQMD parametrization for the exchsve Ai232-/Vj*52o P'^O" 
duction in proton-proton reactions compared to data [56|. The matrix element has been extracted 
from a fit to the exclusive A1232-/V1680 production. 
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FIG. 28. Scheme of a decaying string, ss and uu pairs are created in the color field resulting 
in a hyperon, a kaon and a pion. 




FIG. 29. The Field-Feynman fragmentation function (solid line) is used for newly produced 
particles. A Gaussian fragmentation function (dashed line) is used for leading baryons. 
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FIG. 30. Rapidity spectrum of protons (dotted curve), A's (dashed) and A's (dash-dotted) for 
pp cohisions at 205 GeV/c. Data are taken from [^]. 
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FIG. 31. Rapidity distribution of vr^ and vr in pp collisions at 205 GeV/c. Data are taken 
from m. 
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FIG. 32. da/dxp distribution of protons and A's in pp collisions at 205 GeV/c. Data are taken 
from ||. 
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FIG. 33. Strangeness suppression fs due to a finite transverse string radius R 
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FIG. 34. Pair production rate in a finite color field as a function of the longitudinal field 



extension z p7|| . 
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FIG. 35. Pion 'temperatures' (inv. slope parameter of the pt distribution) extracted from 
e^e~ annihilations in the UrQMD model for different energies are shown. They are compared to 
freeze-out 'temperatures' extracted from a statistical model fit to particle yields [^] in pp, pp and 
e~^e~ reactions. 
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FIG. 36. Mean transverse momentum of vr^ (dotted line) and vr (dashed line) in -pp collisions 
at 205 GeV/c as a function of xp. Data (tt+'s, open squares, and 7r~'s, full squares) are taken from 
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FIG. 37. The mean transverse momentum of vr^'s as a function of number of negatively charg 
hadrons in the reaction p(205 GeV)+p. Note the suppressed zero. Data are taken from |]5^ . 
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FIG. 38. Rapidity distribution of protons and positively charged particles for the reaction 
He+He at y/s = 31 AGeV compared to data p3| . 
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FIG. 39. Dilepton mass spectrum for p+Be at 450 GeV/c. The calculation includes Dalitz 
decays and conversion of vector mesons. Only the curve labeled sum of all contributions (solid 
curve) is folded with the mass resolution of the CERES (full circles) experiment |55| . 
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TABLES 



TABLE I. Baryons and baryon resonances implemented in the UrQMD model. All baryons up 
to 2.25 GeV /(? as well as their antiparticles are included. 



A^(Nucleon) 


A (Delta) 


A (Lamb da) 


S (Sigma) 


SfXi) 


(Omega) 


-^938 


A1232 


A-iiig 


^1192 


"1317 


^1672 


-^1440 


^1600 


■^1405 


^1385 


■-^1530 




A1520 


^1620 


-^1520 


^1660 


"1690 




^ ^ iooo 


^1700 


-^1600 


^1670 


'—'1820 




/V 1 ficn 
^ ^ iboU 


^1900 


-^1670 


^1775 


'—'1950 




/ V 1 C7C 

^ ^ ioTo 


^1905 


-^1690 


^1790 






/V 1 CiQC\ 


^1910 


■^1800 


^1915 






A"i7oo 


A 1920 


A18IO 


S194O 






-^1710 


^1930 


Ai820 


S2O3O 






-^1720 


^1950 


Al830 








-^1900 




A21OO 








A"i990 




A2110 








-^2080 












-^2190 












-^2200 












-^2250 
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TABLE II. Parameters of the CERN-HERA fit |T|l used in UrQMD for the total and elastic 
cross-section above the resonance region {piab > 2 GeV/c). The cross sections are parametrized 
as: <Ttot,ei(p) = A + Bp" + C In^(p) + D ln(p), with the laboratory momentum p in GeV/c and the 
cross-section a in mb. 



a 


A 


B 




C 


D 


n 


pp (total) 


48.0 


0. 




0.522 


-4.51 


0. 


pp (elastic) 


1L9 


26.9 




0.169 


-1.85 


-1.21 


pn (total) 


47.3 


0. 




0.513 


-4.27 


0. 


pp (total) 


38.4 


77.6 




0.26 


-1.2 


-0.64 


pp (elastic) 


10.2 


52.7 




0.125 


-1.28 


-1.16 


7P (total) 


0.147 


0. 




0.0022 


-.017 


0. 


TT+p (total) 


16.4 


19.3 




0.19 


0. 


-0.42 


TT+p (elastic) 


0. 


11.4 




0.079 


0. 


-0.4 


7r~p (total) 


33.0 


14.0 




0.456 


-4.03 


-1.36 


TT^p (elastic) 


1.76 


11.2 




0.043 


0. 


-0.64 


K+p (total) 


18.1 


0. 




0.26 


-1. 


0. 


K^p (elastic) 


5.0 


8.1 




0.16 


-1.3 


-1.8 


K+n (total) 


18.7 


0. 




0.21 


-0.89 


0. 


R-p (total) 


32.1 


0. 




0.66 


-5.6 


0. 


K~p (elastic) 


7.3 


0. 




0.29 


-2.4 


0. 


K'n (total) 


25.2 


0. 




0.38 


-2.9 


0. 


TABLE III. 


Baryon-baryon cross-sections 


in [m^ 


b] from the Additive Quark Model. NN scat- 


tering is explicitly treated, i.e. 


-y/i-dependent, 


etc. 








Bi B2 


N 


A 






n 




N 


40.0 


34.7 




29.3 


24.0 




A 


34.7 


30.0 




25.4 


20.8 






29.3 


25.4 




21.5 


17.6 




n 


24.0 


20.8 




17.6 


14.4 




TABLE IV. 


Meson-baryon cross-sections : 


in [mb] from the Addit: 


Lve Quark Model. 


MB scat- 


tering in the resonance region 


(^ < 1.7 GeV) 


is explicitly treated. 






Ml B2 


N 




A 








TT 


26.6 




23.1 




19.6 


16.0 


K 


21.3 




18.5 




15.6 


12.8 




16.0 




13.9 




11.7 


9.6 
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TABLE V. Meson-Meson cross-sections in [mb] from the Additive Quark Model. MM scatter- 
ing in the resonance region (-^s < 1.7 GeV) is expUcitly treated. 



Ml M2 TT K $ 

TT VL% 142 10?7 

K 14.2 11.4 8.5 

$ 10.7 8.5 6.4 



TABLE VL Mixing angles of meson multiplets according to the flavor SU(3) quark model: these 
parameters assign the pure uu, dd, ss,to the physical particles according to the SU(3) quark model. 
The flavor mixing angles are chosen according to quadratic Gell-Mann-Okubo mass formula p9| . 



For the scalar mesons this formula is not applicable, here an ideal mixing angle (tan(^) = l/\/2) 

is assumed. 

Multiplet degree 

scalar 35 

pseudoscalar -10 

vector 39 

pseudovector 51 

tensor 28 



TABLE VIL Particle multiplicities from the UrQMD per inelastic pp event at 12 GeV/c. Data 
are taken from 



Particle UrQMD Exp. Data 

^ L22 1.44±0.02 

TT- 0.64 0.71±0.02 

0.019 0.019±0.001 

p 1.38 1.27±0.02 

A 0.025 0.037±0.001 



TABLE VIII. Particle multiplicities from the UrQMD per inelastic pp event at 205 GeV/c. 



Data are taken from [52|. 



Particle 



Exp. Data 



UrQMD 



Particle 



Exp. Data 



UrQMD 



TT 

A + So 
p 



2.62±0.06 
3.34±0.24 
0.18±0.05 

0.096±0.01 
1.34±0.15 



2.57 
3.11 
0.16 
0.16 
0.16 
1.32 



K+ 
p 



3.22±0.12 
0.28±0.06 

0.17±0.01 
0.013±0.004 
0.05±0.02 



3.10 
0.26 
0.24 
0.20 
0.037 
0.06 
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TABLE IX. Particle multiplicities from the UrQMD per inelastic pp event at ^/s = 27GeV. 
Data are taken from [58|. 



Particle 


Exp. Data 


UrQMD 


Particle 


Exp. Data 


UrQMD 


7r+ 


4.10± 0.26 


3.79 




3.87± 0.28 


3.72 


vr^ 


3.34± 0.20 


3.16 


K+ 


0.33± 0.023 


0.31 


K- 


0.22± 0.015 


0.22 




0.23± 0.015 


0.26 


V 


0.39± 0.075 


0.36 


Po 


0.385±0.056 


0.50 




0.552±0.129 


0.52 


P~ 


0.355±0.091 


0.41 


UJ 


0.39± 0.026 


0.47 


K*+ 


0.132±0.018 


0.13 




0.088±0.013 


0.080 




0.119±0.023 


0.123 




0.09± 0.017 


0.081 


4> 


0.019±0.002 


0.009 


/2(1270) 


0.092±0.013 


0.119 


P 


1.20±0.119 


1.32 


p 


0.063±0.003 


0.088 


A + SO 


0.125±0.016 


0.19 


A 




0.15 


SO 




0.041 


K + tP 


0.020±0.005 


0.047 


A 




0.038 


TP 




0.009 


S+ 


0.048±0.019 


0.050 


TP 




0.041 






0.0041 






0.0053 


s- 


0.013±0.009 


0.015 


A++ 


0.218±0.016 


0.235 


AO 


0.141±0.019 


0.197 


A++ 


0.013±0.005 


0.016 


AO 


0.034±0.009 


0.026 


S*+ 


0.020±0.004 


0.040 


S*o 




0.071 


T*- 


0.010±0.003 


0.009 
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